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Upper-division student difficulties with the Dirac delta function 
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Department of Physics, University of Colorado, 390 UCB, Boulder, CO 80309 

The Dirac delta function is a standard mathematical tool that appears repeatedly in the under¬ 
graduate physics curriculum in multiple topical areas including electrostatics, and quantum mechan¬ 
ics. While Dirac delta functions are often introduced in order to simplify a problem mathematically, 
students still struggle to manipulate and interpret them. To characterize student difficulties with 
the delta function at the upper-division level, we examined students’ responses to traditional exam 
questions and a standardized conceptual assessment, and conducted think-aloud interviews. Our 
analysis was guided by an analytical framework that focuses on how students activate, construct, ex¬ 
ecute, and reflect on the Dirac delta function in the context of problem solving in physics. Here, we 
focus on student difficulties using the delta function to express charge distributions in the context of 
junior-level electrostatics. Common challenges included: invoking the delta function spontaneously, 
translating a description of a charge distribution into a mathematical expression using delta func¬ 
tions, integrating 3D or non-Cartesian delta function expressions, and recognizing that the delta 
function can have units. We also briefly discuss implications of these difficulties for instruction. 

PACS numbers: 01.40.Fk 


I. INTRODUCTION 

Research into student difficulties at the upper-division 
level is a relatively new but growing area of Physics Ed¬ 
ucation Research [l]. Students in upper-division courses 
are asked to manipulate increasingly sophisticated math¬ 
ematical tools as they grapple with advanced physics con¬ 
tent. To explore this dynamic, some of the literature 
around student difficulties at this level has focused on the 
use of mathematical tools and techniques during physics 
problem-solving . For example, one mathematical tool 
that appears repeatedly throughout the undergraduate 
physics curriculum is the Dirac delta function (hereafter 
referred to as simply, the <5-function). 

The (5-function represents an interesting topic to in¬ 
vestigate student difficulties for several reasons. Despite 
often being discussed as a purely abstract, mathematical 
construct, it is rarely used in purely mathematical con¬ 
texts. In fact, physicists often introduce the (5-function 
into a problem in order to describe or model a concrete 
physical system. Additionally, it is the authors’ experi¬ 
ence that the 5-function is often perceived by experts as 
trivial to manipulate, and is often introduced to simplify 
the mathematics of a problem. However, we still observe 
common and persistent student difficulties using the 15- 
function. Moreover, while the 5-function has appeared in 
literature around student difficulties with Fourier trans¬ 
forms and measurement in quantum mechanics M , we 
are unaware of any existing literature directly targeting 
student difficulties with the 5-function. 

A physics student may be exposed to the 5-function 
at several different points in the undergraduate curricu¬ 
lum. These include its use as a tool to express volume 
mass or charge densities, to describe potential energies or 
wavefunctions in quantum mechanics, in the context of 
Fourier transforms, and, for more advanced situations, in 
the context of Green’s functions. In mathematics courses, 
the 5-function is rigorously defined as a distribution Q, 


and is typically not seen by students until senior or mas¬ 
ters level analysis courses. Such courses are not often 
taken by physics majors. Here, we will focus almost ex¬ 
clusively on the 5-function as a tool to describe volume 
charge densities in the context of upper-division electro¬ 
statics as this is one of the earliest and perhaps simplest 
uses. Despite the general sense of many faculty at the au¬ 
thors’ institution that students will be familiar with the 
5-function by the time they reach the junior-level, we find 
that for some (arguably many) students junior-level elec¬ 
trostatics will be their first exposure to the 5-function. 
Given its many possible uses, we do not claim that the 
research presented here will span the space of all possible 
difficulties with the Dirac 5-function; however, it will pro¬ 
vide a sampling of the kinds of challenges that students 
face when manipulating 5-functions. 

In this paper we utilize an analytical framework Q 
describing the use of mathematical tools in physics prob¬ 
lem solving to structure our investigation and analysis of 
student difficulties with 5-functions (Sec. |TJ)- We then 
present our findings, including common difficulties we 
identified in our student population and a brief discus¬ 
sion of implications for instruction f Sec. Mil) . We end with 
limitations and future work (Sec. II VI) . 


II. METHODS 

Problem solving at the upper-division level is often 
considerably longer and more complex than what is typ¬ 
ical of the introductory level. Making sense of stu¬ 
dents’ work around these upper-division problems can be 
difficult because students often make multiple mistakes 
at different stages of a problem, which then propagate 
through their solutions in unpredictable ways. To help 
manage this complexity, we make use of an analytical 
framework known as ACER (Activation, Construction, 
Execution, Reflection) to scaffold our analysis of student 
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difficulties with the (5-function 00 - 

A. The ACER Framework 

The ACER framework organizes the problem-solving 
process into four general components: activation of 
mathematical tools, construction of mathematical mod¬ 
els, execution of the mathematics, and reflection on the 
results. These components were identified by studying 
expert problem solving 0 and are consistent with both 
a resources 0 and epistemic framing 0 perspective on 
the nature of knowledge. However, the specific details of 
how a given mathematical tool is used in upper-division 
problem solving is often highly dependent on the context 
in which that tool is being used. For this reason, the 
ACER framework was designed to be operationalized for 
specific mathematical tools in specific physics contexts. 
The operationalization process results in a researcher- 
guided outline of key elements in a correct and complete 
solution to a particular problem or set of problems. This 
process will be discussed in greater detail in Sec. Ill Cl 
and additional details about the ACER framework can 
be found in Ref. 0 . 

B. Study Context 

Data for this study were largely collected from the 
first half of a two semester Electricity and Magnetism 
sequence at the University of Colorado Boulder (CU). 
This course, called E&M 1, typically covers electrostatics 
and magnetostatics (i.e., chapters 1-6 of Griffiths 0 )- 
The student population is composed of junior and senior- 
level physics, astrophysics, and engineering physics ma¬ 
jors with a typical class size of 30-60 students. At CU, 
E&M 1 is often taught with varying degrees of interactiv¬ 
ity through the use of research-based teaching practices 
including peer instruction using clickers jl2:] and tutorials 
0 - 

We collected data from three distinct sources for this 
investigation: student solutions to instructor designed 
questions on traditional midterm exams (N=372), re¬ 
sponses to one question from the multiple-response Col¬ 
orado Upper-division Electrostatics (CUE) Diagnostic 
0 (N=146 at CU, N=162 at external institutions), and 
two sets of think-aloud interviews (N=ll). 

Exam data were collected from 7 semesters of the E&M 
1 course taught by 5 different instructors. Of these, four 
were traditional research faculty, and one was a physics 
education researcher. Two of these instructors, including 
the PER faculty member (SJP), taught the course twice 
during data collection. Questions on the exams were de¬ 
veloped solely by the instructor for that semester; how¬ 
ever, in all cases, the question provided a mathematical 
expression for a charge density and asked for a description 
and/or sketch of the charge distribution (e.g., Al-type 
prompt below, see Sec. Ill Cl for the naming convention). 


In one case, the students were also asked to calculate the 
integral of the given (5-function expression. 

Al-type Prompt: 

Sketch the charge distribution: 

p{x,y,z) = cS(x - 1) 

Describe the distribution in words too. 

What are the units of the constant, c ? 

The multiple-response CUE is a research-based, end- 
of-semester conceptual assessment. Only one question 
on the CUE deals with (5-functions, and it is of the same 
general type as the questions asked on the exams (Al- 
type). However, rather than an open-ended prompt, the 
CUE asks a two-part, multiple-choice question. The first 
part provides the student with the mathematical expres¬ 
sion that corresponds to the mass density of two point 
masses (i.e., p(r) = miS(r — rq) + m2i5(r — E 2 )) and asks 
what the integral of this expression is over all space. The 
second part asks what this mass density represents phys¬ 
ically. This question along with the full CUE instrument 
can be accessed at Ref. [15|. CUE data from CU were 
collected from four semesters of the E&M 1 course; two 
of these were courses for which we also have exam data. 
In addition to the CU data, we also collected multiple- 
response CUE data from 9 courses at 7 external insti¬ 
tutions. These institutions range from small liberal arts 
colleges to large research universities. 

Think-aloud interviews were conducted in two sets per¬ 
formed roughly a year apart. Interviewees were paid 
volunteers who had successfully completed E&M 1 one 
to two semesters prior and responded to an email re¬ 
quest for research participants. All interviewees took the 
course during one of the semesters for which we have 
exam data. Final course grades for interview partici¬ 
pants ranged from D to A. Interview protocols for both 
sets of interviews were designed, in part, to probe more 
deeply into difficulties identified in the exam solutions, 
and thus included one or more exam style (i.e., Al-type) 
questions. In addition, the interview protocols also in¬ 
cluded questions designed to target aspects of the ACER 
framework not accessed by the exams and CUE. These 
questions will be discussed in greater detail in the next 
section. 


C. Operationalizing ACER 

The process of operationalizing ACER is discussed in 
detail in Ref. 0. Briefly, operationalization is accom¬ 
plished through a modified form of task analysis 0 0 
in which a content expert works through the problems 
of interest while carefully documenting their steps and 
mapping these steps onto the general components of 
the framework. The resulting outline is then shared 
with other content experts and refined until consensus is 
reached that the key elements of the problem have been 
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accounted for. This expert-guided scheme then serves 
as a preliminary coding structure for analysis of student 
work and, if necessary, is further refined to accommodate 
unanticipated aspects of student problem-solving. 

To guide our data collection and analysis, we oper¬ 
ationalized ACER for problems involving the use of 6- 
functions to express the volume charge (or mass) densi¬ 
ties of 1, 2, and 3D charge distributions. For example, 
the volume charge density of a line charge running par¬ 
allel to the z-axis and passing through the point (1,2,0) 
can be expressed as p(f) = XS(x — 1 )S(y — 2), where 
A is a unitful constant representing the charge per unit 
length. Expressing volume charge densities in this way is 
often necessary when working with the differential forms 
of Maxwell’s Equations and can facilitate working with 
the integral forms of both Coulomb’s Law and the Biot- 
Savart law. 

The elements of the operationalized ACER framework 
for these types of (5-functions problems is described be¬ 
low. Element codes are for labeling purposes only and 
are not meant to suggest a particular order, nor are all 
elements always necessary for every problem. 

Activation of the tool: The first component of the 
framework involves identifying (5-functions as the appro¬ 
priate mathematical tool. We identified two elements in 
the form of cues present in a prompt that are likely to 
activate resources associated with 5-functions. 

Al: The question provides an expression for volume 
charge density in terms of 5-functions 

A2: The question asks for an expression of the volume 
charge density of a charge distribution that includes 
point, line, or surface charges 
We include element Al because, in electrostatics, 5- 
functions are often provided explicitly in the problem 
statement, effectively short-circuiting Activation. 

Construction of the model: Elements in this com¬ 
ponent are involved in mapping the mathematical ex¬ 
pression for the charge density to a verbal or pictorial 
representation of the charge distribution or vice versa. 
Cl: Relate the shape of the charge distribution to the 
coordinate system and number of 5-functions 

C2: Relate the location of the charges with the argu¬ 
ment (s) of the 5-function(s) 

C3: Establish the need for and/or physical meaning of 
the unitful constant in front of the 5-function 
For problems that also require integration of the 5- 
function (e.g., to find total charge from p(r )) there are 
additional elements in construction related to setting up 
this integral. 

C4: Express a differential volume element consistent 
with the geometry of the charge distribution 

C5: Select limits of integration consistent with the dif¬ 
ferential volume element and region of interest 
Alternatively, for sufficiently simple charge distributions, 
one can bypass setting up and computing integrals (el¬ 
ements C4 & C5) by using the physical meaning of the 


unitful constant to simply state the total charge (e.g., 
Q(uniform, spherical shell) = a * 47 tR 2 ). 

Execution of the mathematics: This component of 
the framework deals with elements involved in executing 
the mathematical operations related to the (5-function. 
Because this component deals with actually performing 
mathematical operations, these elements are specific to 
problems requiring integration of the (5-function. 

El: Execute (multivariable) integrals that include one 
or more (potentially multidimensional) (5-functions 
When the results of the integrals in El must be simplified 
for interpretation, Execution would include a second el¬ 
ement relating to algebraic manipulation; however, none 
of the integrals included in this study elicited or required 
significant algebraic manipulation. 

Reflection on the result: The final component in¬ 
cludes elements related to checking and interpreting as¬ 
pects of the solution, including intermediate steps and 
the final result. While many different techniques can be 
used to reflect on a physics problem, the following three 
are particularly common when dealing with (5-functions. 

Rl: Check/determine the units of all relevant quantities 
(e.g., Q, p, the unitful constant) 

R2: Check that the physical meaning of the unitful con¬ 
stant is consistent with its units and the units of all 
other quantities 

R3: Verify that the value of the charge in a region is 
consistent with expectations 
While the first two elements are similar, we consider el¬ 
ement R2 to be a distinct and potentially more sophisti¬ 
cated reflection task in that it is seeking consistency be¬ 
tween the student’s physical interpretation of the unitful 
constant and its units. The distinction between elements 
Rl and R2 was motivated in part by preliminary anal¬ 
ysis of student work, which demonstrated that the link 
between the units of the constant and its physical inter¬ 
pretation may not be explicit for many students (see Sec. 

mm- 

The operationalization of the ACER framework for 5- 
functions revealed that the standard questions typically 
asked on midterm exams (Al-type) do not capture all as¬ 
pects of problem-solving with the (5-function. In particu¬ 
lar, this type of question bypasses Activation at anything 
more than the most surface level assessment of whether 
a student recognize the (5-function. To probe Activation 
more deeply, we began the first set of think-aloud inter¬ 
views with a question that provided a description of the 
charge distribution and asked for a mathematical expres¬ 
sion for the charge density (e.g., A2-type prompt below). 
This A2-type interview prompt also provides a differ¬ 
ent perspective on Construction than the Al-type exam 
prompt by requiring students to generate a multidimen¬ 
sional (5-function expression rather than just interpreting 
it. However, participants in the first interview set often 
failed to activate the 5-function in response to the A2- 
type prompt (see Sec. Mill and thus never moved on to 
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deal with Construction. In order to target the Construc¬ 
tion component more clearly, the second set of interviews 
directly prompted students to use 5-functions to express 
the charge density. 

A2-type Prompt: 

What is the volume charge density of an 
infinitely long, linear charge distribution 
running parallel to the z-axis and passing 
through the point (1,2,0)? Be sure to 
define any new symbols you introduce. 

Sketch this charge distribution. 

Additionally, only one of the seven exam questions 
asked students to integrate an expression containing a 
5-function. To further probe students ability to integrate 
the 5-function (element El), both interview protocols in¬ 
cluded questions that prompted students to calculate the 
total charge within a finite region of space. The second 
set of interviews also ended by asking students to perform 
a set of context-free integrations of various 5-function ex¬ 
pressions (below) in order to more clearly investigate the 
Execution component. These particular integral expres¬ 
sions were each designed to target a specific difficulty we 
anticipated students might have with the procedural as¬ 
pects of integrating the 5-function. While all integrals 
were presented without a physical context, in the case of 
integral d) students were also asked if they could come 
up with a physical situation in which they might set up 
this integral. 

Context-free Integral problems: 

— OO 

a) f 5{x)dx 

OO 

— OO 

b) f xS(x)dx 

OO 

10 

c) f [a5(x — 1) + bS(x + 2)\dx 

o 

d) fff ad(r — r')r 2 sin(6)drd(j)d0 

III. FINDINGS 

This section presents the identification and analysis of 
common student difficulties with the Dirac 5-function or¬ 
ganized by component and element of the operationalized 
ACER framework (Sec. Ill Cl) . 

A. Activation of the tool 

Elements Al and A2 of the framework represent two 
different types of of prompts that can cue students’ to 
activate resources related to the 5-function. In the case 
of A2-type prompts, the student must first recognize that 
the 5-function is the appropriate mathematical tool be¬ 
fore they can correctly answer the question. However, for 


Al-type questions, the 5-function is given as part of the 
prompt, effectively short-circuiting Activation and pro¬ 
viding little information about students ability to recog¬ 
nize when the 5-function is appropriate. 

None of the instructor-written exams included A2-type 
questions. Instead, this element was specifically targeted 
during the first of the two interview sets. When presented 
with the A2-type prompt shown in Sec. Ill Cl 2 of 5 inter¬ 
viewees spontaneously suggested using 5-functions. The 
remaining three participants all expressed confusion at 
being asked to provide a volume charge density for a line 
of charge. Two of these students attempted to recon¬ 
cile this by defining an arbitrary cylindrical volume, V, 
around the line charge and using p = Q/V. This strat¬ 
egy, while incorrect, represents a reasonable attempt to 
make sense of the problem in lieu of the 5-function. Later 
in the interview, all five of these students were presented 
with the 5-functions expression for a linear charge density 
and all but one correctly interpreted this expression as 
describing a line charge. This result suggests that even 
immediately after completing a junior-level electrostat¬ 
ics course, many students may have difficulty recogniz¬ 
ing when the 5-function is the appropriate mathematical 
tool even when they are able to provide a correct physical 
interpretation of it. 

Three participants in the second interview set had just 
completed one semester of upper-division quantum me¬ 
chanics at the time of the interview. To investigate the 
context-dependent nature of student Activation of the 5- 
function, we began these students’ interviews by asking 
for a mathematical expression for the potential of a fi¬ 
nite square well in the limit that the well became very 
narrow and very deep. Two of these three participants 
spontaneously suggested the 5-function as the appropri¬ 
ate tool. Both of these students explicitly focused on the 
use of the words “ very narrow and very deep” just before 
suggesting the 5-function. The third student initially at¬ 
tempted to write the potential as a piecewise function but 
brought up the 5-function when explicitly told there was 
a more compact way to represent the potential without 
using a piecewise notation. Successful activation of the 5- 
function in the quantum case seemed to be linked to the 
high degree of similarity between the visual representa¬ 
tion of a deep, narrow potential well and the commonly 
used graphical representation of 5(x). 


B. Construction of the model 

The Construction component deals with mapping be¬ 
tween the physics and mathematics of a problem. In the 
case of 5-functions, this mapping can be done in two di¬ 
rections: from mathematics to physical description, or 
from physical description to mathematics. The Al-type 
prompts (see Sec. Ill BD used on exams, the CUE, and 
in interviews investigated students’ ability to translate 
a mathematical expression for the charge density into a 
physical description of the charge distribution. As part of 
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this process, students needed to connect the coordinate 
system and number of 5-functions in the given expression 
to the shape of the charge distribution (element Cl). For 
example, the expression p(x,y,z) = cS(x — 1) contains 
one Cartesian (5-function and thus represents an infinite 
plane of charge. Roughly a quarter of students (23%, 
N=87 of 372) had an incorrect shape on exams. The most 
common error was misidentifying volume charge densities 
that included 1 or 2 (5-functions as point charges (62%, 
N=54 of 87). 

On the CUE diagnostic administered at the end of 
the semester, the fraction of students selecting the in¬ 
correct shape was roughly a third (35%, N=51 of 146). 
This trend of selecting the incorrect shape is slightly 
more pronounced in student populations beyond CU who 
have taken the multiple-response CUE (49%, N=79 of 
162, 7 institutions) The most common error for this 
question was misidentifying the given point mass den¬ 
sity (p(r) = miS 3 (r — rq) ) as a solid sphere or spher¬ 
ical shell with radius rq (76%, N=99 of 130 incorrect 
responses, all institutions). This finding is interesting 
because, given students’ tendency to misidentify charge 
densities as point charges on exams, we might expect that 
they would be more successful at identifying point charge 
densities. However, we hypothesize that the prevalence 
of this particular error may have been exacerbated by 
the apparent similarity between the expressions for this 
point mass density and that of a spherical shell (i.e., 
p(r) = rri\6(r — rq)), which are distinguished only by 
the presence of the cube on the 5-function and vectors in 
the argument. 

In interviews, 9 of 11 participants correctly identified 
the shape of one or more charge distributions from the 
mathematical expressions for the charge density. The re¬ 
maining 2 students both sketched the charge distribution 
on 3D Cartesian axes as a very narrow spike originating 
at a point consistent with the argument of the 5-function 
and extending upwards parallel to one of the axes (e.g., 
see Fig.[T]). This ‘spike’ representation was also seen in a 
small number of the exam solutions (5%, N=17 of 311). 
The spike drawn by these students is highly reminiscent 
of the ID graphical representation of S(x) as an infinitely 
tall and thin Gaussian distribution at x = 0 that is com¬ 
monly used when first defining the 5-function. Students 
who draw this ‘spike’ when sketching the charge distribu¬ 
tion may be attempting to apply this ID representation 
of the 5-function to a 3D sketch. The examples shown in 
Fig. [T] have a number of additional interesting features. 
For example, the student in Fig. QJd) places the spike as 
z = 2c rather than z = 2. Additionally, the spikes are 
not always parallel to the same axis despite all having the 
same prompt. These features may warrant further inves¬ 
tigation in future studies; however, as they were observed 
in only a small number of exam solutions and in none of 
the interviews, we will not discuss them in further detail 
here. 

The interviews also included an A2-type prompt (see 
Sec. Ill Cl) to explore students’ ability to translate a phys¬ 


a) 7-j 

C ' 

1 

1 

J2D 

c) 

A 

exi'jU 

y 

Z 




FIG. 1. Examples of the ‘spike’ representation used by 5% of 
exam students and 2 of 11 interview students. In this case, 
students were prompted to sketch the charge distribution de¬ 
scribed by p(x, y, z) = c5(x)5(z — 2). 


ical description of the charge distribution into a mathe¬ 
matical expression for the charge density. This process 
requires students to use the geometry of the charge dis¬ 
tribution to select an appropriate coordinate system and 
number of 5-functions. Of the eight interview partici¬ 
pants who attempted to use 5-functions in response to 
the A2-type prompt, three were able to correctly express 
the line charge density as the product of two ID Carte¬ 
sian 5-functions. Four of the remaining five students in¬ 
stead attempted to use a single 5-function whose argu¬ 
ment was the difference between two vectors, for example 
p oc 5(r — rq) where rq = (1,2, z). These students did 
not explicitly acknowledge this as a 3D 5-function either 
verbally or by cubing the 5-function. Three of these stu¬ 
dents also described the line charge as a continuous sum 
of point charges and integrated their 5-function expres¬ 
sion over all z. This strategy reflects a fundamentally 
reasonable physical interpretation of the situation that 
can be used to construct a correct expression for the line 
charge density; however, none of these students was able 
to fully leverage this interpretation to correctly express 
the charge density. 

Determining the location of the charge distribution (el¬ 
ement C2) was not a significant stumbling block for stu¬ 
dents in this study. None of the interview students and a 
tenth of the exam students (10%, N=39 of 372) drew an 
incorrect position for the distribution. The most com¬ 
mon errors were switching the signs of the coordinates 
(36%, N=14 of 39, e.g., locating the plane described by 
p(x,y,z) = cS(x — 1) at x=-l) or having the wrong ori¬ 
entation of line or plane distributions (38%, N=15 of 
39). However, all relevant exam questions in this study 
dealt with locating charge distributions described by 5- 
functions in Cartesian coordinates, and it is possible that 
student difficulties with element C2 would be more sig¬ 
nificant for other geometries or more abstract notation. 
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The third element in Construction relates to the need 
for a unitful constant in the expression for p(r). For 
Al-type questions, this constant is provided as part of 
the given expression for the charge density, and fully in¬ 
terpreting this expression requires recognizing the phys¬ 
ical significance of this constant. For example, in the 
expression p(x,y,z) = cS(x — 1), the constant, c, rep¬ 
resents the charge per unit area on the surface of the 
plane. Roughly a quarter of the exam students (27%, 
N=70 of 255) presented with an arbitrary constant spon¬ 
taneously commented on its physical meaning and most 
of these (81%, N=57 of 70) had a correct interpretation. 
This fraction should be interpreted as a lower bound as 
additional students may have recognized the constant’s 
physical significance but did not explicitly write it down 
on the exam. The interviews suggest that a students’ in¬ 
terpretation of the constant can be facilitated or impeded 
by their identification of its units. This dynamic will be 
discussed in greater detail in relation to the Reflection 
component fSec. IIII Dl) . 

The A2-type prompts used in interviews, on the other 
hand, do not include the unitful constant but instead 
require that a student recognize the need for this con¬ 
stant independently. Two of eight interviewees did not 
spontaneously include a multiplicative constant in their 
expression for the charge density. Alternatively, three of 
the remaining participants recognized the need for a mul¬ 
tiplicative factor but expressed this factor as p(r) (e.g., 
volume charge density = p(f)5(x — 1 )S(y — 2)). These 
students appeared to be interpreting the p(r) term as 
representing only the magnitude of the charge density 
rather than the magnitude at a specific point in space. 
This also suggests that these students are not treating 
p as the quantity defined by convention to be volume 
charge density and are instead treating it as the quan¬ 
tity conventionally defined as A (i.e., line charge density). 
Ultimately, four of the eight participants included a mul¬ 
tiplicative constant and articulated a correct physical in¬ 
terpretation of that constant before moving on from this 
question. 

The final two elements in Construction (elements 4 & 
5) are not specific to questions involving (5-functions, but 
rather apply to any physics problem involving multivari¬ 
able integration. These elements deal with expressing the 
differential volume element and selecting limit of integra¬ 
tion. We have previously examined student work around 
these two elements in the context of Coulomb’s law inte¬ 
grals [1] and observed a number of difficulties. However, 
only a small number of students (21%, N=12 of 56) strug¬ 
gled with either the Cartesian differential volume element 
or limits of integration in the single semester where the 
exam question asked for total charge in a finite region 
of space. Furthermore, only one of eleven interviewees 
had clear difficulties expressing the differential charge el¬ 
ement in spherical coordinates, and one other selected 
incorrect limits for their Cartesian integral. Thus, set¬ 
ting up integrals for the relatively simple Cartesian and 
spherical geometries used in this study was not a signifi¬ 


cant challenge for our upper-division students. 


C. Execution of the mathematics 

The Execution component of the framework deals with 
the procedural aspects of working with mathematical 
tools in physics. The exams provide limited insight into 
this component as students were asked to actually cal¬ 
culate an integral of a (5-function in only one semester. 
In this case, the students were given the expression for 3 
point charges and asked to calculate / p{f)dr. Roughly 
a quarter of the students (27%, N=15 of 56) made sig¬ 
nificant mathematical errors related to the (5-function 
while executing this integral (element El). The most 
common error (73%, N=ll of 15) amounted to a varia¬ 
tion of equating the integral of the (5-function with the 
integral of the zero point of its vector argument (e.g., 
f S 3 (r — ri)dr = f ridr). Similarly, more than a third of 
students (36%, N=53 of 146 at CU, 41%, N=66 of 162 at 
external institutions) selected an incorrect value for the 
integral of a point mass density in response to the CUE 
question. The most common error (81%, N=43 of 53 at 
CU, 68% N=45 of 66 at external institutions) was equat¬ 
ing the integral of the (5-function with the value (vec¬ 
tor or magnitude) of r at which the argument was zero 
(e.g., J5 3 (r — ri)dr = ri). This issue is different from 
(though potentially related to) the most common issue 
seen in the exam; however, none of the response options 
on the multiple-choice CUE question matched the out¬ 
come of the incorrect setup seen commonly on exams. 
Ultimately, these types of errors indicate that students 
have correctly recognized that the (5-function somehow 
picks out the value r = rq but apply this reasoning in an 
ad hoc fashion in the integration process. 

The first set of interviews investigated Execution in 
the context of calculating the total charge on a uniformly 
charged spherical shell. Two of the five participants were 
not able to complete this calculation and explicitly stated 
this was because they could not recall how to integrate 
the (5-function. When asked for the value of f_ oo S(x)dx 
these participants guessed it would be x or 1/x. Extend¬ 
ing this logic into three dimensions, this response appears 
consistent with some of the common incorrect responses 
on the CUE and exam questions (e.g., f <5 3 (r — ri)dr = rq 
or l/|rq|). 

The second interview set (N=6) targeted the first ele¬ 
ment in Execution differently by asking students to per¬ 
form the context-free integrations shown in Sec. Ill C'1 Two 
students stated that the integral in b) would be equal 
to x without evaluating this expression at x = 0, but 
none of the six participants had difficulty with the in¬ 
tegrals in a) or c). Moreover, the 3D (5-functions used 
on the exam and CUE seemed to evoke different, and 
potentially more fundamental, difficulties than the ID <5- 
functions used in interviews. Three of six interviewees 
also evaluated the r integral in part d) in the following 
way, / (5(r — r')r 2 dr = ^r' 3 , despite correctly executing 
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parts a)-c). Two of these students explicitly stated that 
the effect of the (5-function was only to pick out the value 
r = r'. This result is consistent with that from the ex¬ 
ams and CUE, and suggests that a significant fraction 
of our upper-division students have internalized the idea 
that the (5-function picks out a particular value of the 
variable, but do not recognize the other impacts of the 5- 
functions on the result of the integral, particularly when 
dealing with 3D or non-Cartesian 5-functions. 


D. Reflection on the results 

The Reflection component deals with aspects of the 
problem-solving process related to interpreting and 
checking intermediate steps and the final solution. For 
the questions used in this study, one of the most pow¬ 
erful tools available for checking and interpreting the 
various 5-function expressions is looking at units (ele¬ 
ments R1 and R2). In particular, for Al-type problems, 
looking at the units of the given constant (e.g., the c 
in p(x,y,z) = c5(x — 1)) can facilitate interpretation 
of that constant’s physical meaning. Five of the seven 
exam prompts explicitly asked students to comment on 
the units of the given constant and two-thirds of the stu¬ 
dents (70%, N=178 of 255) responded with the correct 
units (element Rl). 

Beyond just commenting on the units of the given con¬ 
stant, we would like our students to recognize the physi¬ 
cal interpretation of this constant (element C3), and en¬ 
sure that the units and physical interpretation are con¬ 
sistent (element R2). However, it was often difficult to 
assess if students did this on the exams, in part because 
only a quarter of our students (N=70 of 255) explicitly 
commented on the physical interpretation of the unit¬ 
ful constant (element C3, see Sec. IIII Bl) . Most of these 
students (83%, N=58 of 70) also provided units for the 
constant that were consistent with their physical inter¬ 
pretation of its meaning. However, a third of the students 
(32%, N=82 of 255) gave units that were inconsistent 
with the shape they identified regardless of whether they 
commented on the constant’s physical meaning. Some of 
these students had the correct shape but incorrect units 
(40%, N=33 of 82), and others had the incorrect shape 
but correct units (30%, N=25 of 82). For an expert, 
the units of the constant, its physical meaning, and the 
shape of the charge distribution are tightly linked; how¬ 
ever, these results suggest that this relationship may still 
be developing for the students. 

Interviews offer additional insight into the connection 
between the units and physical interpretation of the unit¬ 
ful constant. When prompted to comment on units, 9 of 
11 participants explicitly argued (incorrectly) that the 
5-function is unitless and thus, regardless of the geome¬ 
try of the charge distribution, the units of the constant 
must be C/m 3 . This argument was often justified by the 
statement that the 5-function was ‘just a mathematical 
thing’ and thus did not have units. The belief that the 


5-function is unitless was surprisingly persistent, and in 
some cases interfered with students’ physical interpreta¬ 
tion of the unitful constant. For example, four of these 
nine students had articulated a correct physical argument 
for the units of the constant prior to arguing that the 5- 
function was unitless. Two of these students abandoned 
their previous physical argument in favor of a unitless 5- 
function, and the other two were unable to reconcile their 
intuition about the constant’s units with their belief that 
the 5-function should be unitless. Ultimately, only 2 of 
the 9 students were able to independently convince them¬ 
selves that the 5-function had units of 1/m (in this case). 
The other 7 students were only convinced when specifi¬ 
cally prompted by the interviewer to consider the units 
of dx in the expression f S(x)dx = 1. 

The third element in Reflection deals with using the 
expression for the charge density to calculate the total 
charge within a finite region in order to ensure this value 
is consistent with expectations. While this is something 
that we saw experts do spontaneously, almost no students 
(N=3 of 255) explicitly executed such a check on the ex¬ 
ams when not prompted to calculate total charge. Sim¬ 
ilarly, only one of the interview students spontaneously 
attempted to calculate total charge in a region to gain 
confidence in his answer. Interview students were also 
explicitly prompted to calculate total charge and to con¬ 
sider the units of this expression (element Rl) later in the 
interview. This process helped 4 of the 9 students who 
argued that the 5-function was unitless to realize there 
was a problem with the units of their expression. These 
findings underscore the potential value of these reflective 
practices (element Rl-3) both in terms of catching errors 
and facilitating interpretation, but also suggest that our 
students rarely engaged in reflection spontaneously. 


E. Implications for Instruction 

This study was not designed to investigate the impacts 
of different instructional strategies or curricular materi¬ 
als on the prevalence or persistence of students’ difficul¬ 
ties with the Dirac 5-function; however, our findings do 
suggest several implications for teaching or using the 5- 
function at the upper-division level. For example, in¬ 
structors should be aware that students are unlikely to 
encounter the 5-function in their required mathematics 
courses, and that it may or may not be covered in a math 
methods course run within a Physics department. Thus, 
the assumption that all students in a junior-level course 
will be aware of the 5-function and its properties may not 
be justified. Students’ struggles with the procedural as¬ 
pects of integrating more complex 5-function expressions 
may be one manifestation of this lack of sufficient prior 
experience with the 5-function. These students may ben¬ 
efit from additional opportunities to practice integrating 
both ID and 3D delta functions and in multiple coordi¬ 
nate systems. 

Additionally, canonical 5-functions questions rarely, if 




ever, require a student to consider when the (5-function is 
an appropriate tool. Questions that describe a charge dis¬ 
tribution and ask for an expression for the charge density 
can provide a baseline assessment of student ability to ac¬ 
tivate the (5-function when not prompted explicitly. This 
type of question also addresses another finding, that con¬ 
structing a mathematical expression for the charge den¬ 
sity is a distinct and potentially more challenging task for 
our students than interpreting that same expression. As 
the former task is arguably the more authentic, students 
may benefit from additional opportunities to construct 
various 5-function expressions in multiple coordinate sys¬ 
tems. 

The belief that the 5-function is unitless was a sur¬ 
prisingly prevalent and persistent idea. This belief may 
be exacerbated by presenting the 5-function as a purely 
abstract mathematical construct. Moreover, the idea of 
a unitless 5-function can interfere with students’ inter¬ 
pretation of the unitful constant. To facilitate student 
reflection on problems involving the 5-function, specific 
emphasis should be placed not only on the fact that the 
5-function can have units, but also on how to determine 
them based on the argument of the 5-function. 

IV. CONCLUSIONS 

This paper contributes to the limited body of research 
on student difficulties with the Dirac 5-function by pre¬ 
senting an application of the ACER framework to guide 
analysis of student problem-solving with the 5-function in 
the context of expressing charge densities mathematically 
in junior-level electrostatics. The ACER framework pro¬ 
vided an organizing structure for our analysis that helped 
us identify nodes in students’ work where key difficul¬ 
ties appeared. It also informed the development of inter¬ 
view protocols that targeted aspects of student problem 
solving not accessed by traditional exams, particularly 
around activating the 5-function as the appropriate tool 
and executing integrals of various 5-function expressions. 

Our upper-division students encountered a number of 
issues when using/interpreting the 5-function. These dif¬ 
ficulties represent a subset of students’ difficulties with 
the 5-function and may not include issues that might 
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